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SECOND GENUS CROSSED ORBITS! 


By DANIEL BUCHANAN? 


Abstract 


In a former paper, first genus periodic orbits were obtained for the motion 
of the electrons in the “‘crossed orbit” model of the normal helium atom as 
proposed by Kemble and Bohr. The present paper treats the second genus 
orbits, i.e., orbits which are in the vicinity of the first genus orbits and which 
re-enter after the lapse of many first genus periods. 

The methods employed involve the solutions of differential equations with 
periodic coefficients. Notation is introduced which greatly expedites the 
obtaining of these solutions. 


1. Introduction 


The problem considered in this paper is to determine second genus periodic 
orbits for the motion of two infinitesimal bodies that repel each other and are 
attracted by a third finite body. The forces of attraction and repulsion vary 
according to the Newtonian law of the inverse square. 

As the three bodies in this problem are somewhat analogous to the normal 
helium atom, the finite body will be referred to as the nucleus and the infini- 
tesimal bodies as the electrons. 


2. First Genus Orbits 


One periodic orbit is that in which the electrons move in a circle with the 
nucleus as centre, the electrons remaining diametrically opposite. This plane 
of motion will be called the equatorial plane. Angular distance on either side 
of it will be called latitude, while longitude will have the obvious meaning. 
In a previous paper (2) the author obtained two types of three-dimensional 
periodic oscillations in the vicinity of this circular orbit. In the first case the 
longitudes differ by 180°, whereas the latitudes are equal in magnitude but 
opposite in sign. In the second case the longitudes differ by 180°, but the 
latitudes are equal in magnitude and also in sign. 


The results that were obtained in the paper just cited are similar to those 
that were first proposed by Kemble and later independently by Bohr for the 
“Crossed Orbit’’ model of the normal helium atom (6, Chap. VII). 


The present paper treats the second genus orbits for Case I of the author’s 
previous paper cited. 


1 Manuscript received September 13, 1941. 


Contribution from the Faculty of Arts and Science, The University of British Columbia, 
Vancouver, B.C. Read at the May, 1940, meeting of the Royal Society of Canada. 
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3. The Differential Equations 


Choose a system of rectangular axes with the nucleus as origin and the 
equatorial plane as the xy-plane. Let the ratio of repulsion to attraction be 
denoted by k? and let the units of time and of distance be so chosen that the 
gravitational constant of attraction is unity. If the co-ordinates of the 
electrons are x1, 1, 21, and x2, ye, 22, the force function is 


jet, 2. 2 


T1 12 T2 
mn = (xi + yi + 2), Yo = (xo + yz + 23)", 
[(x1 — x2)? + (yr — ye)? + (21 — 22)*)%, 


and the differential equations of motion are 


ll 


Ti2 














ax; ou dy; o0U as; _ aU (1) 
dt? 7 OX; dt? OY; dt? 02; 
j= 1,2. 
In Case I of the paper cited above 
“a = — i, VWF —- VW, BSE KRM. (2) 
Then 7, = re = 3r2. It is then necessary to consider the motion of only 
one electron, say 1, 1, 21, and the differential equations of its motion are 
Cn wm, En hn, Ga _ me, (3) 
dt? ri dt? ni dt? ni 
where w = 1 — 4k? and 0 < uw <1. Transforming to cylindrical co- 


ordinates r, 0, z and ae, use of the integral of areas n@ =C, 
dt dt 


a constant, we obtain the equations 





a 

dz r (72 + 22)8/2" 

o i (4) 
de OF + 


4. The First Genus Orbits 


A particular solution of Equations (4) is 
r= c/p, z = 0, 


being the case in which the electrons move in a circle in the equatorial plane 
and remain diametrically opposite. To obtain the differential equations of 


displacement put 
< (1 + €p) “ (et) 
r=— Ep), z = —(€), 
Me . be 


\ 
| 
a (5) 
b= fom 5 +O)", 
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where p, ¢ are new dependent variables, Tr is the new independent variable, 
and 6 and € are parameters. Thus 


p+ (1+ 6)p = R = (1+ 4) | (3p +30) — 6€(p*? + pf?) 


+ (0p + 15p'%? — Po + --- |, (6) 
E+ (+ dg =Z=(1 + 8) | 3epo — e6ptg —F 59) + & (lop*s — 2 ps) 
4. et ( = 15944 +2 pegs - 25) ss we ret, |: 


where the dots denote derivation with respect to t. The right member R 
is even in ¢ and Z is odd in ¢. In the coefficients of the various powers of € 
it is found that in R the parity of p is opposite that of €, and in Z it is the same 
as that of €«. From these properties it is seen that p is odd in € while [ and 6 
are both even in e. If the initial values 

60) = $0) =0, 0) =1, (7) 


are taken, the periodic solutions of Equations (6) are 


; (3 + cos 2T)e — a (51 + 12 cos2r + cos 47) +--, 


p = 
¢ = sint + (O)e + (O)et+ ( Deo + —--, (8) 
6 = 32+ (O)e* + ( De® + --. 


When these results are substituted in Equations (5) the equations for the 
periodic orbits of the first genus are obtained. 


SECOND GENUS ORBITS 
5. Definition 


Given a set of differential equations 


dx; nad . sae aearite 
“dt nen X s(x; , €; t), 7 = i, ,n, 


where the X; are analytic in the arguments x;, € and do not contain ¢ 
explicitly. Suppose they admit the periodic solutions 


x: = 6,(e, t) 
having the period JT. These solutions are said to be of the first genus. If 
we substitute a ee sig Sak ikea Pe 
in the differential equations then solutions for the y;, say 

yi = Wi (€0,A;2), 
can be found having the period 

NT(1 + a power series in A), 
N being an integer. The solutions 

x; = O:(€, t) + Vi(eo, A; A 
are said to be of the second genus (5, Vol. 3, Chap. XXVIII). The second 
genus orbits approach those of the first genus as the parameter A approaches zero. 
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6. The Differential Equations 

Let € = € (1 + A), where € occurs explicitly in Equations (6) but not 
where it enters implicitly through p, ¢, and 6. Let po, £0, do represent the 
values in Equations (8) when € is replaced by €). Further, let 

p = pot P, fS=f+¢ (9) 

and let (1 + 6) be replaced by (1 + 6) (1 + Y) where #, g are dependent 
variables and ¥ is a parameter that vanishes with A. The above substitutions 
will be made in Equations (6) and for the time being the subscript o on e, 6, 
p, ¢ will be dropped. Hence Equations (6) become 

p+p+(1+ d(1t+y)(o +p) = (1+ 8) + [3d + AY {(o +)? 

+A +9") — 6e(1 + N*E(0 + P+ (0+ PF +9} +---], 

F+qt+ Hd tyF+a = 1+ H+ y[3ed + (+ B) (10) 

(§ +g) — 321 + A/{2Ap+ PPG +q —AE+O}4+---] 

where the terms independent of ~, g, y, \ cancel off. When the solutions 
for p, ¢ are substituted in the above we obtain 


b+ Pup + Prg = P,, 


q = Qup + Qiog = Or ; (11) 


where 
3 3 
a 3 5 e(1 — 3cos 27) + 35 (17 + 20 cos 27 + 27 cos 47) + ---, 
Pe = Qu = — 3esinT + 5.e(ts sin T — sin 3r)+ --, 
Qxe = 1- z e(1 — cos 2T) + 3e4(1 — cos 2T) + --. 


The values for P; , Q, will be determined later. 


7. The Equations of Variation 


If the right members of Equations (11) are neglected, 


p + Pup + Prog = 0, q + Qup + Qeug = 0, (12) 
and these are called the equations of variation. Their generating solutions are 
1 1 
a ee ol Se ae 
€p a8 (3 + cos 27) 64° (51 + 12 cos 27 + cos 4r) + (13) 


ef esint +( )e®+ -—---. 
To obtain the solutions of Equations (12) the usual substitutions 


p a en, q = ey, (14) 


are made, where @ is a constant and u and v are functions of rt. Equations (12) 


then become 
(D? + 2iaD — ao? + Py)u + Pov = 0, \ 
0 


Qiu + (D? + 2iaD — a? + Qi2)v ) 


’ 
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The existence proof, which is omitted here, shows that a, as a power series 
in €, and certain initial values can be determined so that u and v will be 
periodic with the period 27. 

Now Pi: and Qi are odd in € while P;,, Qi2 and @ are even in €. Then u 
and v are of opposite parity in €. Since Solutions (14) are later multiplied 
by arbitrary constants we may for one set of solutions select u as even and 
v as odd or u as odd and v as even. The former selection will be made and 
accordingly we assume 


soa + a? oe + a? e+ - - <, 
v= wmMety7% B+ ---, (16) 
a=1+me@+ameée+—--. } 


Let Equations (16) be substituted in Equations (15) and let the resulting 
equations be denoted 1s (15’). We then equate to zero the coefficients of the 
various powers of € in (15’) and we thus obtain sets of differential equations 
that define the various u“” and v°*, It will now be shown that the various 
Q@; and certain constants of integration can be determined so that each 
u?? and v%+» will be periodic. 

The question of the initial values for u and v will be deferred until after 
the first step of the integration. 


Step 1. Terms in (15’) Independent of « 
There is but one differential equation at each step. At this step it is 

(D? + 2iD)u™ = 0, 

u™ = Ao + Boe, 
There are two arbitrary constants A» and By. If we assign the initial value 
u(0) = 1, the solution would be uw = Ay + (1 — Ao)e'7. The integration 
was carried out with this initial value and the u’s and v's were found to contain 
exponentials having the highest multiples of 7 exceeding by two the power 
of € associated with the u’s and v’s. This leads to a somewhat complicated 
form for the final solutions. We therefore choose the initial value “(0) = 0 
from which it follows that each u°°(0) = 0. The desired solution at the 
first step is then u = Ao, where Ao remains arbitrary. 


and its solution is 


Step 2. Terms in € 
The differential equation is 
(D? + 2iD)v) = (3 sin 7)u, 
and the solution expressed in terms of exponentials is 
® = C+ Dew + oer + 3e-7), 


where Ay, C,, D; are arbitrary. 


Step 3. Terms in & 
The differential equation is 


(D? + 2iD)u® = U® = (20, - 3 a : cos 2) w+ (3 sin T)v™. 


2 
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Non-periodic terms will arise from the constant terms and terms in e~*7 in 


U™, These terms are 
U® = 2A oe ~ (0)e-7. 


Hence either Ao or @ must be zero. To put Ayo = 0 would lead to the trivial 
solution, hence we must put a, = 0. The integration then yields 


u = A, — : Affe + ¢™) + [Cie* + 3(C, — Die* + Dye-*7], (17) 
where Ay, Az, Ci, D; are arbitrary. 


Step 4. Terms in & 
The differential equation is 
(D? + UDyp® = VO = — : (11.sin r —.sin 37)u + (Gaine)a® 


— 3 (1 — cos 27)v™, 
Now V® contains the terms 
y@ =3 (< se ay) is (sc, -4 Ds) enti, (18) 
and these terms must be annulled if v is to be periodic. Hence C, = D, = 0 


5 ; ; ‘ ‘ 
and u® = A, — 3 A,(e* + e'7), The integration for v® will then be 


8 
where Ao, Az, C3, Dz are arbitrary. The constants C3, Ds; remain arbitrary 
‘through Step 5 but at Step 6 they enter V™ the same as C,, D, occur in 
Equation (18). Hence €; = D3 = 0 and these terms will be ignored in v® 


above. 


vo = Cy + Dse2it + ae (ei? + 3e-7) — 1A o (3e7 — 3e-i7 — de-sir), 


Step 5. Terms in e4 
The differential equation is 
(D? + 21D)u® = U® = 2au™ — ; (1 — cos 27T)u 


a < (17 + 20 cos 27 + 27 cos 47) u™ (19) 


+ onary < : (tunr —sadee™. 


So far as the constants and terms in e-?'7 enter U™ we have 
99 3 45 ; 

To on is ree Se aes A Cees See —2iT 

U Ao( 20 ) (54: ig 40) € , 





32 
Hence 
99 15 
a = — Gj Ao, A, = = Ao, 
and the desired solution for u“ is 
335 s BF 93 5 165 
oO ae en meee” dE Pack OF... et 
. A Ao (S32 ¢ a ~ jee” ae }. 


where Ao, Ag are arbitrary. 
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The General Step 
At the 7th step, U“ contains the terms 


U® = Ao(2as + constant) — {5 Ac + Ao (constant) } ce, 


and there are no other terms where a and A, enter. Hence as and A, can 
be uniquely determined. The arbitrary constants arising through a v-equation 
enter the next v-equation as in Equation (18) and must be annulled. An 
induction to the general term will readily show that the a’s and the constants 
of integration can be uniquely determined as above. Since Ao remains 
arbitrary we may take Ao = 1, as we later multiply u by an arbitrary constant. 

One set of solutions of the equations of variation has thus been found, 
viz., 


p = ery, , q= e770, , (20) 
where 
uy = 1 — + (cos 27) + ---, 
oo J 
= 2 2 [AY” cos 2kr + 1.BY sin 2kr]e, 
j=0 k=0 
% = €(2costT—Ztsint)+( )@+—-—, 
© j 
= 2 2 [Cf’* cos (2k + 1) r + 7 D&t” sin (2k + 1)7]e%*?, 
j3=0 k=0 


A second set of solutions can be obtained by changing the sign of 7 in 

Equations (20) thus 

pb =e" um, q = &" (— iw), (21) 
where a 

uw = um, vg = Vy. 

The remaining two sets of the solutions of the equations of variation can 
be obtained, as shown by Poincairé (5, Vol. 1, Chap. IV), by partial differenti- 
ation of the generating solutions Nos. (13) with respect to the two arbitrary 
constants that they contain, viz., the initial time ¢) and the scale factor e. 


Differentiating Equations (13) first with respect to tp we obtain 


ee ae ee a a ee 
dto a (1+ 6)* ar a (i+ 6)**’ 
where 
uz = — 5 (sin are + az (6 sin 27 + sin 4r)e* + ---, 
2 ee ee ee 
Ito ec (1+ 6)* ar ce (1+ 6)* *’ 
where 


v3 = (cosT/et+( )?@+-—-. 
Since these solutions are later multiplied by arbitrary constants, we may drop 
the constant multipliers of us; and v; and take 
p = Me, gq = % (22) 


as the solutions. 








18 CANADIAN JOURNAL OF RESEARCH. VOL. 20, SEC. A. 


Differentiating Equations (13) with respect to € we obtain 


i (2¢22) d(ep) a7 96 


de OT «06 Oe’ 


(22) def) O7 96 


q* de de 


OT 06 de’ 


where the parentheses about the partial derivatives denote partial differentia- 
tion only in so far as € enters p and ¢ explicitly. Performing the differentia- 
tions and multiplying the results by € so as to make the power of € the same 
as that of the highest multiplicity of 7 in its coefficient, we obtain 


Pp = um+ Kru, g = 4+ K7v;, 
ly = ; (3 + cos 27) — ig (51 + 12 cos 27 + cos 4r)e* + --, 


(sint)e +( JDP +H+—-—-, 


= oe 9 
K = mS ae ee 


” 


V4 


I 


The criterion that the four sets of solutions that have been obtained con- 
stitute a fundamental set is that the determinant formed from these solutions 
and their first derivatives must be different from zero. ‘This determinant is 


uy Us Us us + Kus 
Uy + 1QU; Uz — iQ, U3 us + K(tTuzs + U3) 
A=|. ; (23) 
1V — We V3 vy + KTv3 
| iv. — an, — iW. — A’. 3 vg + K(Tv3 + 3) 


Since the equations of variation do not contain 7 explicitly, this determinant 
is a constant (4, §18). Its value may be obtained with the minimum of com- 
putation by putting 7 = 0. We thus obtain 
A.= — 2ie [1 + & (a power series in e*)], 

and this is different from zero for € not zero but sufficiently small numerically. 
Hence the four solutions constitute a fundamental set, and the most general 
solutions of the equations of variation are 

bp = Kieu + Ke eu, + Kgus + Ka (us + Krus), 

eS 1(K, ev, —_ Ke e7*7y2) ad K3v3 a K, (v4 + Krvs) > 
where K,, K:, K;, Ky, are arbitrary constants. 


(24) 


8. Notation 


As we shall have frequent occasions in the sequel to use power series in odd 
or even powers of € with sums of sines or cosines of odd or even multiples 
of 7 in the coefficients, we shall adopt a notation* that will represent these 
series and denote their properties. 

The letters C and S will denote the power series referred to above, C 
representing series with cosines in the coefficients and S with sines. Their 


*This notation was first used by the author in his paper ‘‘Pertodic orbits of the second genus 
near the straight line equilibrium points in the problem of three bodies” (1). 





BUCHANAN: SECOND GENUS CROSSED ORBITS 19 


properties will be indicated by superscripts consisting usually of two paren- 
theses. In the first parentheses will appear a digit followed by the letter o 
or e. The digit will represent the power of the lowest term in ¢€, and the 
parity of this digit will be the same as that of the parity of € in the series, 
zero being taken as even. The letters o and e will signify that the multiples 
of 7 in the sines or cosines are odd or even respectively. 

The highest multiple of 7 in the coefficient of €& is usually 7 but occasionally 
it exceeds 7. When it is the same as j no second parenthesis in the super- 
script will appear. When it exceeds j a second parenthesis will be added 
containing a digit denoting the excess of the multiple over the power. Thus 


S 


« j 
Ce = > (as? cos 2k) &, 
j=0 E=0 


© 1 
Uae « = 
i=0 k= 


s 


it 
z= (agit? cos 2kr) et!, 
j k=0 
© j 
SOOD = DY Z (bFA,sin (2k + 1)r) e. 
j=0 k=0 

The symbols C and S with superscripts but without subscripts will denote 
types of series only. When a particular series of a type is to be noted, sub- 
scripts will be used. 

Expressed in the above notation, the generating solutions and the w’s 


and v’s of the equations of variation are represented as follows: 


p = Cow, C= Siow, 
uw = C%) + 7§%9, “a, = i, 

1 = Ce) 4 7500, v = v1, 

us = SZ, v3 = Ci, 
te = C2, m% = §2., 


9. The Highest-Multiple-Power Property 


Certain pairs of series like us, Us have the property that the coefficient of 
the sine of the highest multiple of 7 in the coefficient of any power of € in 
one series is the same as, or differs only in sign from, the coefficient of the 
cosine of the highest multiple of 7 in the coefficient of the same power of € 
in the other series. This property will be called the Highest-Multiple-Power 
property and will be abbreviated to H.-M.-P. The word same or opposite 
will follow according as the coefficients in question are of the same or opposite 
signs, respectively. Thus 

us, us, H.-M.-P. opposite 
vz, v4, H.-M.-P. same. 


10. Construction of the Periodic Solutions 


We proceed to show that Equations (10) can be integrated as power series 
in \ and that y and the constants of integration can be uniquely determined 
so that p and g will be periodic with the period 7, a multiple of 27, and will 
satisfy certain initial conditions. 
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Let us assume 
> = p> pid’, pa pa qi’, 4S >) vit’, (25) 
and let these substitutions be made in Equations (10). The resulting equa- 
tions will be cited as (10’). On equating the coefficients of the various powers 
of \ in Equations (10’) we obtain sets of differential equations from which 
the p; and the q; will be found. It will be shown that the y; and the arbitrary 
constants arising at each step can be so chosen as to satisfy the periodicity and 
certain initial conditions. 
As the first genus orbits are symmetric, that is, 6(0) = ¢(0) = 0, we 
shall confine our attention to the symmetric second genus orbits. Hence 


we put p(0) = g(0) = 0 and when these conditions are imposed upon 
Equations (25) we have 
i) = q;(0) = 0, Gj = 1, 2, 3, = —). 


It will be found that, under the above initial conditions, one constant of 
integration will remain arbitrary at each step of the solutions for p; and q;. 
We may therefore impose an additional condition and we shall choose p(0) = 1. 
Hence #:(0) = 1, $;,(0) = 0,7 >1. There is no loss of generality here as p 
carries the factor \ in Equations (25). The parameter \ then becomes the 
scale factor of the second genus orbits. 


Step 1. Coefficients of \ in Equations (10') 
On equating the coefficients of \ in Equations (10’) we obtain 
pr + Pupi + Pog = PY = 7,02 4+ COM, 


qi --- Outi “bh Qi0q1 = ow = 7150.0 — S20 | (26) 
The complementary functions of these equations are 
Pi = RM! M Uy + RP eT, + RSPus + RY (us + Kus), (27) 
Qa = (RPE, — key.) + ksPvs + RY (ve + Kvs), 
where k{”, — — -, ki? are the arbitrary constants. By using the method of 
the variation of parameters to find the particular integrals we have 
RMeary, + RMeary, + Rus + RY (us + Kus) = 0, 


Mela (4, + iam) + Re (tt, — icrus) 

+ RP tis + RY {aug + K(rus + us)} = P, (28) 

RP ein, — RPE “Tv. + RsPust RY (vs + Krv3) = 0, 

het (G5, — a) — be (td, + arve) 

+ PIs + RY {ds + K(ris + v3)} = Q. 

The determinant of &{”, — — -, R{? is the same as in Equations (24) and 
therefore does not vanish for € ~ 0 but small numerically. Equations (28) 
are therefore solvable, the solutions being 


Aki = — e@T(MyP™ + MuQ”), 
Ak” = eT (MoP® + MnQ®), - (29) 
AR® = — (MxP”® + MuQ), 


ARY = (MeP® + MuQ®), 
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where Mu(j = 1, 2, 3, 4; k = 2, 4) are the co-factors of the elements in 
the determinant A, j referring to the column and k& to the row. The forms 
of these co-factors in terms of the notation adopted, together with their 
properties, are given in Table I. 

TABLE I 


THE CO-FACTORS AND THEIR PROPERTIES 




















Co-factors Type of series Functions H.-M.-P. 
properties 
My = Mn (C2) + 1.812%) 
| 
Mu = Mu 2E2(Cig) + 7Sy 4+) 
Me | i(C§° + KrS9°) | 
CS? SE? Opposite 
Me 453° 
Mss 4( Si -e KrCg) 
3 a, Co | Same 
Mus | Co” 





On making the necessary substitutions in the first two of Equations (29), 
integrating and substituting the results in Equations (27) we obtain particular 
integrals for ~: , g: which are respectively of the same form as the right mem- 
bers of Equations (26). Non-periodic terms, however, will arise if @ is an 
even integer. Those values of € that will make a@ an even integer must 
therefore be excluded. 

The integration of the last two of Equations (29) proceeds differently from 
the preceding integration. When the series involved are expressed in the 
notation employed, 


pgp = ind? + eds? + cy?” a Cio + Kr(yS$9 + S10), 


a (1,€)(3) 1,€)(3 
RY? = —7,S2Z°° — SZ — y,S0@ — Suo@, 


where 
C32? ,, Ss have H.-M.-P. opposite, 


CY, Si have H.-M.-P. opposite, 
and where the d{” here and in the sequel denote power series in €? with constant 
terms different from zero. 
When the integrations are performed and the results for kj” and kj” 
are substituted in Equations (27), the non-periodic terms arising from the 
integration of 71S;!°° + Sj°® in Rk, R cancel off and we obtain 


A = 710% 4 CBOM 4 KMys + (xe -f. : yd + cd rus, 
se (30) 


1520 a S20) + K3"03 a (xxp a ! vid; a cd‘? ) TV3. 


qi 
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bNRo 
ht 


Because of the H.-M.-P. properties listed, the numbers in the second 
parentheses of the superscripts of C and S are 1 instead of 3. 


From the periodicity conditions we must put 
KK,» + : yd? + eds” = 0, (31) 


and the complete solutions at the first step become 

pi = KPe*u, + KMe yy + Kuz + KOuy + W1C1OD 4 CAOM, 

qa = i(KPe*v, — Ke-*7y») + Kv3 + KPg + 71, S09 + S200, 
where K{”, — — -, Kj” are the constants of integration. The part of the 
particular integrals in Equations (30) obtained from the integration of RD 
and k;” have been absorbed by the corresponding particular integrals arising 
from Rk) and R2. 

The period of the above solutions must be a multiple of 277/a@ and also of 
27. Hence a must be rational, i.e., a = Ni/N, where N, and N are relatively 
prime integers. The period will then be 


T = N, (7) = Nx). (33) 


From the initial conditions £,(0) = 0, or gi(0) = 0, and p,(0) = 1, it 
follows that 


(32) 


K® = K® = 1, K = 0. 
The solutions, therefore, which satisfy the initial and periodicity conditions 
y J 
are iar lfiet 
pi = eu, + eu + KP ug + 1C1® 4 Cao, 


34 
qa = i(e*"v, Pe ey») of. K{v4 + 150) + S20) | ( ) 


where Kj” and ¥; must satisfy Equation (31). 


Step 2. Terms in * 

It will be necessary to consider the terms in \? where another relation 
between K{” and y; will arise which together with Equation (31) will 
uniquely determine these constants. 

The differential equations at this step are 

be + Pute + Peg = P, 

Qe + Quire + Qn@ = Q°, 

and the complementary functions are the same as Equations (27) with 

?, --— —-, kf as the arbitrary constants. We shall require first only 

the terms in e*““’ in the right members. So far as these terms are concerned, 
P® = e* [eK (C8 4 i829) + (CO +4 180%) + CO 4 7520] 


(35) 


+ e"| conjugate ., (36) 
QO? = ie [eK (CO +759) + ¥4(C8 4 7S9%) + C4 4 (S00) 
— ie conjugate ¥ 


In the above and in the sequel the word ‘‘conjugate’’ in the parentheses 
designates the conjugate of the function in the parenthesis immediately 
preceding. 
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On varying the parameters k;”’, k? we obtain 
Ak® — e*[MypP® + MuQ), 
Ak? = e'@ [MoeP® + MuQ]. 
When the substitutions are made for the M’s, P’s, and Q’s in Equations (37), 


constant terms will arise which must be annulled if the solutions at this step 
are to be periodic. This gives the second relation in Kj” and y: 


(37) 


KP + 5 yd? +a” = 0, (38) 
which along with Equation (31) will give 
KY =1d?, ys = ede. (39) 
On substituting Equations (39) in Equations (34) the solutions at the first 
step take the form Pr = eu, + em, + COM 
nie ea’ a (40) 
Q = ev, — ey, + SCO, 


When we make use of Equations (40) instead of Equations (34), the right 
members P”, Q@ take the form 


P® = eetiar(C0O 4 4509) +4 ¢€c-ia (conjugate) 
ee'27(C + 75%) + ee" (conjugate) 
C1.0M 
Q® = eetiar(GCI 4 S00) +4 ¢€e-tiar (conjugate) 

+ e@"(iC% + S$) + e-@"(conjugate) 

+ Seow, 
Varying the parameters k}, kf we have equations similar to the last two 
Equations (29), if we replace the superscript (1)with (2). These two equations 
will give a relation between K? and 2 similar to Equation (31). When the 
integrations are completed and we impose the further condition p.(0) = 0 
we obtain 


- 
+ 


9 
bo = 2 (et@7 UP? + e-iier U}?) + C10, 
j=l 


(eliar V2 — e-iier YD) 4 Seow 
1 


q = 1 


W brs 


J 


where U;P is similar to u,, Vj? to v1, ; UjP , UP are conjugates and so also 


are Vi, V#?. 

The integration for the higher terms in » and g may be carried on as far 
as desired. One relation K{ and y; will be obtained at the step where KY? 
arises and a second relation will be found at the next succeeding step. The 
general form will therefore be 
Pn =_ > (e#o7 Uj? + e—tiar U)) + C10 | 


j=l 


n 
Qn = 4 2 (ella Yim _ e~tiar Vm) + S20 D | 
j=l 
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and the final form for 7 and zis 
Cc © 
oe [it+te(pt+ 2 pad*)I, 
c © 
2=—e > Gnd"). 
: c+ 7 ) 


The period in ¢ of the first genus is 


3 


and of the second 
Tee MRE 2 ed"). 


11. The Convergence of the Solutions 


Only the formal construction of the solutions has been made and it now 
remains to consider their convergence. The usual method of establishing 
the convergence is by an existence proof based upon Poincaré’s extension to 
Cauchy’s theorem (5, Vol. 1, pp. 58-63). Such existence proofs are long and 
actually more difficult than the formal construction. An alternative method 
is to make use of MacMillan’s theorem (3) where he showed that if the 
constants of integration can be determined so as to make the solutions of 
differential equations similar to those arising here formally periodic, then such 
solutions will converge for all finite values of the time, provided the scale 
factor, such as X, is sufficiently small numerically. 


References 


. BucHANAN, D. Proc. Roy. Soc. (London) A, 114 : 490-516. 1927. 

. BUCHANAN, D. Rend. Circolo Mat. Palermo, 60: 1-21. 1931. 

. MacMititan, W. D. Trans. Am. Math. Soc. 13 (2) : 146-158. 1913. 

. Mouton, F. R. Periodic orbits. Carnegie Institution, Washington, D.C. 1920. 

. Porncaré, H. Les méthodes nouvelles de la mécanique céleste. Gauthiers-Villars, Paris. 
1892-1899. 

. VAN VLECK, J. H. Bull. Nat. Reseafch Council. Vol. 10, Pt. 4, No. 54. 1926. 





